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Abstract 

We use Monte Carlo simulations to study a dynamically triangu- 
lated disk with Ising spins on the vertices and a boundary magnetic 
field. For the case of zero magnetic field we show that the model 
possesses three phases. For one of these the boundary length grows 
linearly with disk area, while the other two phases are characterized 
by a boundary whose size is on the order of the cut-off. A line of con- 
tinuous magnetic transitions separates the two small boundary phases. 
We determine the critical exponents of the continuous magnetic phase 
transition and relate them to predictions from continuum 2-d quantum 
gravity. This line of continuous transitions appears to terminate on a 
line of discontinuous phase transitions dividing the small boundary 
phases from the large boundary phase. We examine the scaling of bulk 
magnetization and boundary magnetization as a function of boundary 
magnetic field in the vicinity of this tricritical point. 

1 Introduction 



The dynamical triangulations approach to (Euclidean) 2-d quantum grav- 
ity ^[21 El arose from attempts to apply the Feynman path integral prescrip- 
tion for quantizing field theories to gravity. A straightforward application of 
the prescription would involve integrating over all 2-d manifolds of a given 
topology. For generic matter coupled systems this integration cannot be 
done. The dynamical triangulation method replaces the integral with a sum 
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over discretized manifolds composed of identical equilateral triangles con- 
nected along their edges. In addition to including several exactly solvable 
models, the approach is also amenable to computer simulation using tech- 
niques from statistical physics. This is the approach we have taken in this 
paper. 

Most work based on 2-d dynamical triangulations so far has been re- 
stricted to closed topologies where Liouville theory provides many predic- 
tions. The use of a topology with a boundary allows us to probe less well 
understood areas — in this case the Ising model coupled to a boundary 
magnetic field (H) and 2-d quantum gravity. This has been studied ana- 
lytically [H ij . The boundary magnetic field can be seen as interpolat- 
ing between the only two conformally invariant boundary conditions: fixed 
(H = iboo) and free {H = 0). Away from these special boundary conditions, 
the model possesses the interesting property of being conformally invariant 
in the bulk, but not at the boundary. We describe our model in more detail 
in Section [21 

In principle, the model has four couplings. In the gravitational sector, 
there are couplings to the disk area (number of triangles, A) and to the 
boundary length (/). In the matter sector, there are the Ising coupling (J) 
and the boundary magnetic field. We describe our methods for simulating 
the model in Section EJ They were initially tested against a strong coupling 
expansion of the partition function. This along with measurements of the 
disk amplitude for zero Ising coupling are discussed in Sectional Both agree 
with predictions. 

For fixed disk area and zero boundary magnetic field, we have mapped 
the resultant 2-d phase diagram. We show evidence for three phases and 
associated critical lines, and use finite size scaling to extract estimates for 
critical exponents. These lines of transitions appear to meet at a unique 
tricritical point at which the system exhibits long range spin correlations and 
a non-trivial scaling of boundary length with disk area. This is presented in 
Section El 

As explained in Section ^ we have tracked this tricritical point as a 
function of boundary magnetic field. We show our results concerning the 
scaling behavior of the model with non-zero boundary magnetic field and 
test some of the conjectures contained in [T. Finally, we describe possible 
future work. 
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2 The Model 



Our model consists of a 2-d triangulation composed of equilateral triangles 
connected along their edges with Ising (7j spins on the vertices. We allow 
only combinatorial triangulations; that is we exclude degenerate triangula- 
tions — no two triangles have the same three vertices, no triangle includes 
the same vertex twice and each triangle may have only one neighbor per 
edge. We characterize a configuration by its area, A, boundary length, /, 
and Ising spins, cTj. The action for such a configuration is 

S = — J (TjCTj — H (7fe + fiA + ul (1) 

(ij) ked 

where {ij) are nearest neighbor vertices, d denotes the boundary of the 
disk, J is the Ising coupling, and /i and ly are the cosmological constant 
and boundary cosmological constant respectively. The action contains no 
explicit gravitational component because, in two dimensions, the Gauss- 
Bonnet theorem implies that the Einstein-Hilbert action is a topological 
invariant 0- The grand canonical partition function is 

A,l tAl {crj 

where X^i^; denotes a sum over all allowed triangulations of area A and 
boundary length /. 

3 Simulation Methods 

A disk can be created from a sphere by simply cutting out a hole, and 
this observation forms the basis of our simulation method. If we consider 
a triangulated sphere with a single marked vertex and imagine cutting out 
this vertex and all the triangles containing it, we arrive at a triangulated 
disk. Furthermore every triangulated disk can be reached from an associ- 
ated triangulated sphere. This allows us to sample the disk ensemble by 
sampling from a sphere ensemble using the usual procedures |H] — the only 
modification being that we identify one vertex as special, never allowing our 
local geometry changing moves to delete it. Additionally we do not endow 
it with an Ising spin. The special vertex along with the ring of triangles 
which contain it, form a cap which converts the disk into a sphere. Fig. ^ 
shows an illustration. 
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Figure 1: Making a disk from a sphere. 



Our sampling of triangulations and Ising states is done with a Metropo- 
Us procedure [Hj. Several types of moves are used in the simulation. For 
changes to the geometry, the moves are link flips, vertex insertions and ver- 
tex deletions. For changes to the Ising matter, we use either a local spin 
flip, or a Wolff cluster algorithm pL(]^0 




3.1 Geometry Moves 

To change the geometry, we choose a triangle at random, then choose at 
random one of three types of geometry changing moves: link flips, vertex 
insertions, and vertex deletions. For vertex insertions, this completely spec- 
ifies the move. For vertex deletions, we need to choose which of the chosen 
triangle's three vertices to delete, while for link flips, we need to choose which 
of its three links to flip. Only vertices other than the special marked vertex 
and with exactly three surrounding triangles are allowed to be deleted. No- 
tice that when a new vertex is added we must also assign a new Ising spin. 
Moves are disallowed if they would lead to a degenerate triangulation. 

These moves are ergodic in the space of triangulated spheres, and since 
every disk can be converted into a sphere, they are also ergodic in the space 
of triangulated disks. In addition to ergodicity, we require detailed balance. 
Let f be a triangulation, o the set of spins on t unchanged in the move and b 
the set of spins on t which potentially change in the move. The requirement 
of detailed balance is 



^ Later revisions of our code which use versions of these algorithms updated for the 
Potts model have been used for some of this work. 




e-S(i'^'b)T{t,a,b^t',a,b') 



(3) 
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where S{t, a, b) is the action and T(t, a,b ^ t' , a, h') is the transition proba- 
bihty. Each of our geometry changing updates requires the transition prob- 
abihties satisfy this criteria. We have used two types of geometry move in 
our simulations — a simple Metropolis move and a more efficient heat bath 
style update — which are described below. 

3.1.1 Simple Metropolis Algorithm 

In this method, any spins which will be changed by the move are chosen at 
random from a uniform probability distribution in advance. In that case we 
have h' so 

r(t, a,h^ t', a, h') = Pp{t, a,b^ t' , a, b')Pa{t, a,b ^ t' , a, b') (4) 

where Pp is the probability of proposing the move and Pa is the probability 
of accepting the move. This leads to the following update probabilities: We 
accept a proposed vertex insertion with probability 

Pl{t, a, b - t', a, b') = ^min{e-^(*>'''')+^(*"^'^ 1} (5) 

where A{t) is the number of triangles in t. A proposed vertex deletion is 
accepted with probability 

P^{t\ a, b' ^ t, a, b) = lmin{e+^(*''"'^')-^(*'"'^\ 1} (6) 
while link flips are accepted with probability 

P/(i, a, 6 ^ t', a, b') = min{e-^(*''"'^')+^(*'"'''\ 1}. (7) 

3.1.2 Heat Bath Algorithm 

In the heat bath variation, insertions or deletions are accepted or rejected 
without regard to the state of the vertex's spin. In the case of insertions, 
the vertex's spin is chosen from a heat bath distribution after the move has 
been accepted. The probabilities in this scheme are: 

Pait', a, b' ^t,a)= ^ ^ ^_J^ ^ (8) 

and 

P-(t.a^f,«,6') = ^ ,^^_,„^^_,,. (9) 
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where AS-^ is the change in the action if the involved spin was or will be 
chosen to be up, and ASi if it was or will be chosen to be spin down. The 
spin of inserted vertices are chosen from: 

-S{t',a,b') 

3.2 Spin Updates via Local Spin Flips 

The procedure for updating a single Ising spin is to choose one at random 
and flip it subject to a Metropolis test. The probability of accepting a flip 
is 

P^{a -a) = min{e-(^(-'^)-^('^)\l}. (11) 



3.3 Spin Updates via the Wolff Algorithm 

The Wolff cluster algorithm |10j is an efficient way of updating Ising (or 
Potts) spins. Specifically, it has been shown to be very efficient at combating 
the effects of critical slowing near phase transitions. The algorithm begins 
by selecting a vertex at random. The spin of this vertex is assigned as the 
first spin of a cluster. The cluster is grown by examining all the vertex's 
neighbors and adding to the cluster any with the same spin as the initial 
vertex with probability 1 — e"^"^. The cluster is then grown further by 
following the same procedure with all of the newly added vertices. The 
procedure is iterated until no new vertices are added. Finally, all the spins 
in the cluster are reassigned to a (single) random value. The Wolff cluster 
algorithm satisfies detailed balance for the Ising model: 

P ^ ! ! ! 1 = (12) 

P^{t,a,b' ^t,a,b) e-^^(*''^''') ^ ' 

where Pp is the probability of proposing a move (growing a given cluster) 
via the Wolff update. Since, unlike the Metropolis algorithm, all moves are 
accepted this is the same as the probability of performing the move. The 
cluster algorithm relies on the fact that the action contains only nearest 
neighbor bond couplings. At non-zero boundary field this is not the case 
and we have used a generalization of the algorithm which handles non-Ising 
like terms. 

Consider an Ising model with such non-Ising like additions to the action: 
S{t, a, b) = Si{t, a, b) + So{t, a, b) (13) 
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where Sj is the Ising action and So is the new term. Let Pp{t, a,b ^ t, a, b') 
be the probabihty of proposing a transition from (t, a, b) to (t, a, b') and 
Pa(t,a,b t,a,b') be the probabihty of accepting the proposal. Now we 
require 

Ppjt, a,b^t, a, b') Pgjt, a,b^t, a, b') ^ e-^^(*'°'^') e-^o{t,a,b') 

Pp{t, a, b' t, a, b) Pa{t, a, b' t, a, b) e-^Ht,«,b) e-3oit,a,b) ■ ^ > 

If we ignore the non-Ising interaction and propose transitions as the 
Wolff algorithm prescribes we will have 

Pp{t, a,b^t, a, b') = P^{t, a,b ^ t, a, b'). (15) 

This in conjunction with 1121 and 1141 gives. 

Pa{t,a,b^t,a,b') ^ e-^o(*'"'''') 

Pa{t,a,b' ^t,a,b) e-^o{t,a,h)- ^ > 

We satisfy this with the choice 

Pa(i,a,6^t,a,6') =min{e-(^o(i,a,6')-5o(t,a,b))^l}. (17) 

Thus we modify the usual cluster update so as to subject the entire move to 
a Metropolis step dependent only on the non-Ising like terms in the action. 

For present purposes, the non-Ising addition to the action is the inter- 
action of the boundary spins with the boundary magnetic field: 

So = -HY,<^i. (18) 

3.4 Critical Tuning, Fixing Area and Boundary Length 

We expect the disk amplitude for the fixed area and boundary ensemble to 
be of the form [TTl IH E] 



2 



Z{1, A) = Te''^' AV=^e~ A . (19) 

In the above, (Xc and Vc are the critical values of the cosmological constants. 

Given the simulation methods we have described, both area A and 
boundary length / will fluctuate. We can not prevent these fluctuations 
without sacrificing ergodicity. However, if we modify the action by adding 
a geometry fixing part 

5/(A„ A^, h, Ao) = Xu{l- lo? + \M - ^0? (20) 
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we can force the simulation to remain in the vicinity of some target area 
Aq and boundary length Iq. Notice that these modification terms vanish for 
disks with these target values and we can hope to recover the results of the 
fixed area/boundary length ensemble with A = Aq and / = /q by merely 
sampling such target disks from our ensemble. 

With the above addition, the simulation will visit states with weight 

W{l,A) = e"'^-''(-^'"^^''"'^o)ZM (21) 

2 

= ia^('^c-u)l-X^(l-lo)^ J^b^{^J.^-^J.)A-X^(A-Ao)'^ ^-^ 

If we assume that the average values of / and A will occur approximately 
where they maximize W, we have 



= (22) 



(23) 



which we can solve to get 

'~{Af {A) 



and 

^c = -^ + ^^ + 2A.((0-/o)-^. (25) 

These equations allow us, as a by-product, to read off estimates for the 
critical couplings fic and fc by measuring the mean area (^4) and boundary 
length (/) . Of course to do this, in general, we need the values of the 
exponents a and b — numbers that are not generally known (see Section ^2] 
for a counter example). Notice, however, that the values of a and b are 
irrelevant in the thermodynamic limit (^O)^o — > oo) so long as Iq <C Aq. 
This condition will be met close to the tricritical point where / ~ \/A. 



4 Disk Amplitude and Strong Coupling 

We have tested our simulation methods in two ways. In the first, which 
is akin to the small coupling expansion, we measure the frequency of visits 
to small triangulations and compare the results to known values. In the 
second, we measure pure gravity disk amplitudes and compare estimates of 
the exponents for area and boundary length with known values. 
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A, I 


count 


visits 


1,3 


1.00 


1.00 ± 0.00 


2,4 


1.50 


1.50 ± 0.01 


3, 3 


1.00 


0.99 ± 0.01 


3, 5 


3.00 


2.99 ± 0.01 


4,4 


3.75 


3.75 ±0.03 


4,6 


7.00 


6.96 ±0.10 


5,3 


3.00 


2.99 ±0.01 


5,5 


12.60 


12.60 ±0.10 


5,7 


18.00 


17.93 ±0.20 



Table 1: Comparison of the number of distinct triangulations for a given 
area and boundary to the frequency of visitations by the simulation with no 
Ising spins. 

4.1 Strong Coupling 

We may recast the partition function in the following form 

Z = (wie-^ ± wae"^^ + (Jae"^^ . . .) (26) 

where to a represents the weight associated with triangulations of area A with 
all other variables summed over. For large n, only the smallest triangulations 
survive and the coefficients oja may be calculated explicitly. This is referred 
to as the strong coupling expansion. Computationally, we can achieve the 
same result by explicitly limiting the simulation to triangulations less than 
a certain area. Unlike the above, this procedure does not require large /u. 

We hand counted the number of distinct triangulations for disks with 
areas up to five. The counting is weighted by inverse symmetry factors, the 
calculations of which were the most difficult part of the counting. For the 
case of = u = and no Ising spins, the number of times the simulation 
visits triangulations of a given area and boundary length should be propor- 
tional to this count. In Tab. ^ we compare our hand calculated counts to the 
observed number of visitations by the simulation. The data are normalized 
so that the values in the first row are one. The agreement was excellent and 
served as an important confirmation of the correctness of our code. 

We extended this calculation to triangles with Ising spins and possible 
boundary fields. The results, shown in Tabs. Inland |31 were also excellent. 
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A,l 


count 


visits 


1,3 


1.00 


1.00 ±0.00 


2,4 


3.03 


3.03 ± 0.04 


3,3 


2.04 


2.03 ± 0.04 


3,5 


12.27 


12.24 ±0.30 


4,4 


15.44 


15.43 ± 0.30 


4,6 


57.89 


57.78 ± 1.00 


5,3 


12.45 


12.45 ±0.30 


5,5 


104.94 


104.73 ± 3.00 


5,7 


301.05 


300.61 ± 5.00 



Table 2: Comparison of the number of distinct triangulations for a given 
area and boundary to the frequency of visitations by the simulation with 
Ising spins included with J = 0.1 and H = 0. 



A,l 


count 


visits 


1,3 


1.00 


1.00 ±0.00 


2,4 


5.29 


5.30 ±0.20 


3,3 


2.07 


2.09 ±0.10 


3,5 


37.38 


37.44 ± 2.00 


4,4 


27.17 


27.31 ± 1.00 


4,6 


308.57 


309.93 ± 20.00 


5,3 


12.76 


12.81 ± 1.00 


5,5 


320.80 


321.67 ±20.00 


5,7 


2808.60 


2820.27 ± 100.00 



Table 3: Comparison of the number of distinct triangulations for a given 
area and boundary to the frequency of visitations by the simulation with 
Ising spins, a Boundary Magnetic Field with J = 0.1 and H = 1. 
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4.2 Disk Amplitude 



If we eliminate the Ising spins and tune the cosmological constants to their 
critical values, we obtain the disk amplitude which is predicted to take the 
form dU Ha E] 

Z{A,l) = l^A'>e-^. (27) 

The exponents a and b are predicted to be — | and — | respectively. We have 
attempted to measure these exponents from our simulations. The exponent 
b can be found by measuring the fraction of bulk vertices which are 3-fold 
coordinated, while a can be obtained by fitting the so-called 'baby universe' 
distribution. 



4.2.1 Measuring the exponent b 

Let tl be a triangulation os a disk of area A — 2 and boundary length I. Let 
t2 and ^3 be disks of area A and boundary length /. Detailed balance tells 
us 

W{ti)P{h ^ t2) = W{t2)P{t2 ^ h) (28) 

where W{ti) is the weight of state ti and P{ti — > t2) is the probability of a 
transition from ti to t2- From this we get 

W(h) Pih -^t2) = Y. W{t2)P(t2 ^ tl). (29) 

We can use this to write the partition function as 

Z{A-2,l) = ^Wih) (30) 



t2 



(31) 



Now consider the ratio 

Since we have eliminated the Ising spins, the action is for all states and the 
transition probabilities are 1 for all possible moves. That being the case, the 
sums over transition probabilities just count the number of possible moves. 
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The number of possible insertions is equal to the number of triangles into 
which a new vertex can be inserted: 



Y,P(ti^h) = A-2, (33) 

ts 

while the number of possible deletions is equal to the number of bulk vertices 
with exactly three surrounding triangles, N^,: 

J2P{t2^ti) = N3. (34) 

So, 



Z{A-2,l) ^ I N3 
Z{A,l) \a-2/^(,^) 

Given the expected form of Z, 

N3 \ [A 



^ ~ 2 / V A 



(35) 



(36) 



■■■(:i^>,.,,, - --Mt) 

We generated data with areas between 50 and 1600 to test this and found 
it was statistically consistent with h = —2.5. Because the above partition 
function is only correct for large volumes, the best results are obtained by 
excluding disks smaller than some cutoff from the fits. Of course, excluding 
too much data will harm the fits — however our data shows a robust plateau 
close to the predicted value. This can be seen in Fig. |2l 

4.2.2 Baby Universes and the exponent a 

When two sections of the disk share a boundary of only one link, we say that 
the disk has split into two minimal neck baby universes. Using the form of 
the partition function mentioned above, the probability of a disk of area A 
and boundary / generating a baby universe of area Ai and boundary li is 



li{l + 2-h)Y {Ai{A-A 



I ) \ A 
^ ll {l + 2-hf l\ 
^^P^ ^2 [A-AiY A^^ 



(39) 
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Figure 2: 6 vs. smallest included volume. 



A, I 


a 


a 


Xd.o.i. 


500, 32 


-0.451 


0.0008 


246.2 


1000, 44 


-0.487 


0.001 


41.0 


2000, 64 


-0.498 


0.001 


13.4 



Table 4: a which gives best fit to baby universe histogram. 

We generated histograms of baby universe observations with various fixed 
disk areas and boundary lengths. To do this we found all such minimal 
necks on various equilibrium configurations and accumulated 2-d histograms 
labelled by the areas and boundary lengths of the smaller baby universes. 
Using the analytical value of b = —2.5, we found the value of a which gave 
the best fit of the above probability distribution to the data. These are 
shown in Tab 0] The result for A = 2000 is within 2a of the predicted value 
a = — ^ and we believe that the discrepancies for A = 1000 and A = 500 
are likely finite size effects. 
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5 Phase Diagram with H = 



For the work described in this section, both cosmological constants were 
set to their critical values. Areas were kept fixed, but Aj^ was set to zero, 
allowing the boundary length to fluctuate. We have identified phases based 
on two observables, magnetization and boundary length. To locate phase 
transitions we look for peaks in the corresponding susceptibilities. We define 
the magnetization of a given disk to be 



-4 



(40) 



where N is the number of vertices in the disk, and the sum is over all vertices. 
We use angle brackets to indicate averages over sampled disks, so (M) refers 
to an average of M over disks generated by the simulation. The magnetic 
susceptibility is 

d (M) 

XM = = {NMM) - {NM) {M) . (41) 

When we refer to boundary susceptibility, we mean 

XL = '-^ = {L-)-{L?. (42) 

By these criteria we have found three phases — a paramagnetic phase 
with a small boundary (I), a ferromagnetic phase with a small boundary 
(II) and a paramagnetic phase with a large boundary (III). The locus of 
susceptibility peaks for A = 1000 is shown in Fig. |21 and a schematic phase 
diagram is shown in Fig. |3 

We used finite size scaling to determine the order of the observed tran- 
sitions. To make clear possible finite size effects, we looked at exponents for 
successive pairs of volumes. The scaling exponent for some quantity x with 
values xi at area Ai and X2 at A2 is given by 

\n{x2/xi) 

HM/A^y ^ ' 

For sufficiently large boundary cosmological constant, the boundary 
length is minimal and we expect behavior similar to that of a marked sphere. 
In that case, there should be a continuous magnetic phase transition at 
Jc = — ^ Intanh (^ In ^) ~ 0.216273. The observed magnetization and 
susceptibility in Fig.s El and IHl show this. Continuum calculations for 2-d 
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X Peaks in Boundaiy Susceptibility 
O Peaks in Magnetic Susceptibility 
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0.4 
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^ure 3: Observed peaks in susceptibilities for A = 1000, H 
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Figure 4: Schematic phase diagram at H = 0. 
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Figure 5: Magnetization, 1/ = 5, H = 0. 



0-0 1000 triangles 
B-B 2000 Triangles 
0—0 4000 Triangles 




Spin Coupling 



Figure 6: Magnetic susceptibility, i' = 5, H = 0. 
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A 

Areas 


ln(x2/a;i)/ 111(^2/^1) 


1000 - 2000 


0.70 ±0.01 


2000 - 4000 


0.70 ±0.01 


4000 - 8000 


0.69 ±0.02 


8000 - 16000 


0.68 ±0.02 


16000 - 32000 


0.67 ±0.04 



Table 5: Scaling exponents of x = (^NM'^'^ at transition between phases I 
and II. 



Areas 


ln(x2/xi)/ln(^2Mi) 


1000 - 2000 


0.832 ± 0.006 


2000 - 4000 


0.765 ± 0.008 


4000 - 8000 


0.68 ±0.01 


8000 - 16000 


0.60 ±0.02 


16000 - 32000 


0.51 ±0.04 



Table 6: Scaling exponents of x = (^NM'^'j for u = —5 and J = 0.75, H = 0. 

quantum gravity coupled to c = | matter predict (^NA'P'j ~ As Tab.l^l 
shows values of this exponent extracted from data at succeeding pairs of 
areas. The agreement is quite good for the larger areas. 

For sufficiently small v, the boundary length is maximal and, as shown 
in Fig. [7| the geometry is like that of a branched polymer. In this essentially 
1-d configuration we expect that the system will not magnetize. However, 
from Fig.s|niand|Hl it appears that the line between Phases I and II continues 
into small v. Notice, though, that the susceptibility scaling exponents on 
this continuation, shown in Tab. 101 decrease as the disk area is increased. 
This indicates that the continuation is a finite size artifact. To confirm this 
we examined how the mean magnetization varies with area at u = —5 and 
J = 0.75 which would be within this small-boundary, ferromagnetic phase 
if it existed. We fit the data to the form cq + ci * A^'^ with the constraints 
<= Co, ci <= 1. As shown in Fig.|Hl extrapolation based on this fit implies 
that the magnetization will approach zero as the area is sent to infinity. To 
the extent that the fit does not agree with the data, the fit seems to be 
underestimating the rate of decrease so the predication appears safe. Our 
conclusion is that for small there is only one phase (III) characterized by 
a large boundary and unmagnetized spins. 
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Figure 7: A triangulation, J = 5, = —5, H = 0. 




Spin Coupling 



Figure 8: Magnetic susceptibility, i' = —5, H = 0. 
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40000 



Area 

Figure 9: Magnetization vs. area, J = 0.75, v = —5, H = 0. 



Areas 


ln{x2/xi)/\n{A2/A,) 


1000 - 2000 


1.38 ±0.01 


2000 - 4000 


1.07 ±0.02 


4000 - 8000 


0.97 ±0.02 


8000 - 16000 


0.99 ±0.002 


16000 - 32000 


0.96 ±0.02 



Table 7: Scaling exponent of x = — (/)^ for J = 1.0, v = 1.3 and H = 0. 

Fig. El and Tab. [3 show that along the transition from small (I, II) to 
large (III) boundary phases, the mean boundary length appears to jump 
discontinuously as u is varied and the boundary length susceptibility exhibits 
a critical exponent of approximately one. This is strong evidence that the 
transition is first order. 

6 Scaling Behavior with a Boundary Field 

At the tricritical point, we performed simulations to test the predictions of 
regarding the scaling of the bulk magnetization, the boundary magnetization 
and their susceptibilities with the boundary magnetic field. For this part 
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of the investigation, the Ising couphng was set to its continuum critical 
value. To reduce the number of parameters we needed to keep track of, the 
boundary lengths were fixed at \fA. While this does involve tuning v to its 
critical value, the boundary-length- fixing Gaussian in the action makes its 
exact value unimportant. Had we not done this, the approximate nature of 
the boundary coupling tuning procedure would have made comparing data 
from different runs more difficult. 

The existence of a boundary magnetic field breaks the symmetry between 
up and down spins, so we no longer need to take the modulus of the spin 
sums. We can now define magnetizations 



Ma 



^E'^^ (44) 



N -I 



Ma = yE^- (45) 



The susceptibilities are defined in the usual manner. Expressions for both 
boundary and bulk magnetizations are derived in [1]: 

IMS (e^^-l)(3 + (2 + ^/7)e^^) 
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and 



LiA 3 



(47) 



where 



L = 



a{H)l 



(48) 



and 



a{H) 




(49) 



However, at least some of these results are dependent on the exact na- 
ture of the discretization procedure, i.e. the nature of the lattice ensemble. 
Our ensemble of so-called combinatorial triangulations is a subset of the set 
of triangulations employed by these authors. In concrete terms we might 
expect that the function a{H) relevant to a description of the combinato- 
rial ensemble might differ substantially from the function given above. In 
spite of this it is reasonable to expect qualitatively similar behavior in the 
two cases. We have therefore looked for a boundary magnetization which 
increases with H and is independent of A and /, as well as a bulk magne- 
tization which decreases with H for fixed /. The boundary magnetization 
plot shown in Fig. II II indeed illustrates this predicted behavior. In contrast, 
Fig. 1121 shows the bulk magnetization increases with H. 

Caroll, Ortiz and Taylor speculate that their predicted decrease of bulk 
magnetization with H is due to an increase in the average spin's distance 
from the boundary when the boundary field is increased. To test this, we 
measured geodesic distances from all vertices on the disks to the boundary 
as follows. All boundary vertices are labelled as having a geodesic distance 
from the boundary of zero. All vertices connected to the boundary vertices 
and not already labelled are labelled as having a geodesic distance from the 
boundary of one. We continue labelling successively deeper vertices until 
there are none left. At that point we compute the average distance from 
the boundary. We observed behavior opposite that predicted by Caroll, 
Ortiz and Taylor: the mean distance to the boundary, illustrated in Fig. 1131 
decreases with increasing magnetic field. This behavior is then consistent 
with the behavior of the bulk magnetization. 

An additional prediction is that the boundary magnetic susceptibility 
should be finite, Xd = "^^3^ 1 ™ the limit of ^ and A — > 00. This 
differs markedly from the flat-space behavior of infinite susceptibility jl4j 
which exhibits a logarithmic singularity. To compare this result with the 
numerical simulations we must be careful. The matrix model result is a 
prediction for the behavior of the boundary susceptibility on an infinite 



21 




Figure 11: Boundary magnetization vs. boundary magnetic field at the 
tricritical point. 
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Figure 12: Bulk magnetization vs. boundary magnetic field at the tricritical 
point. 
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Figure 13: Mean geodesic distance from boundary vs. boundary magnetic 
field at the tricritical point. 



area disk as the boundary magnetic field is tuned to zero. Our simulations 
are carried out at finite area and so we should not take H all the way 
to zero to do the comparison (the magnetic field effectively fails to break 
the symmetry of the Ising action for fields H < l/l). Thus we are led to 
examine the dependence of the susceptibility on the disk area for small but 
non-zero magnetic fields. In Tab. |H1 we show finite size scaling exponents 
for the boundary magnetic susceptibility at three small boundary magnetic 
field strengths. While the exponents at H = 0.0216 are positive, which 
would indicate a divergence, this field strength is below the value needed to 
break the spin-up, spin-down symmetry for the smallest disks of area 1000 
(■y=^ ^ 0.0316). For the larger values of H, we find negative exponents 
which imply a finite value of the susceptibility for small magnetic field in the 
infinite area limit. This is in qualitative agreement with the matrix model 
predictions and quite different from the flat space result. 

We have looked at the behavior of three geodesic correlation functions 
which measure the average properties of the spins as a function of distance 
to the boundary. The first of these, ni(r, N), measures the mean number of 
vertices encountered at distance r from the boundary. The second, ^^(r; A'^), 
measures the net (integrated) magnetization and the third, M{r;N), the 
mean magnetization at distance r from the boundary. They are defined 
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H 


A1,A2 


lose Y9 / Yi 

logh/h 


0.0216 


1000,2000 


0.0863 




2000,4000 


0.0891 




4000,8000 


0.0651 


0.108 


1000,2000 


-0.0190 




2000,4000 


-0.0218 




4000,8000 


-0.0197 


0.216 


1000,2000 


-0.0490 




2000,4000 


-0.0361 




4000,8000 


-0.0636 



Table 8: Critical exponent of boundary magnetic susceptibility at the tri- 
critical point. 

explicitly as follows 

ni{r;N)= (E^'^A.r) ~ iV^-^/'^«Fi(r/iVV'^«) (50) 
na(.r;N)= {E^<^iSD,,r) N^-"-^'"" F^{r /N^'"") (51) 
M{r-N)= (E.f^) ~iV-V/(r/iV'/'-) (52) 

where we have also introduced finite size scaling forms for these three func- 
tions analogous to the functions used for geodesic correlation functions de- 
fined on the sphere in Requiring that the area {N ~ A) dependence in 
M(r; A^) disappear for small r implies that this correlation function takes 
on a simple power law form for small x 

gM{x) ^ x-"""- (53) 

so that 

M{r]N) r~'^H^. (54) 

Figs. 11411151 and ll61 show our finite size scaling data for a representative value 
oiH = 0.433 and yield estimates for l/du = 0.27 ±0.05 and A = 0.1 ±0.02. 

The value for the fractal dimension, dn ~ 4, is consistent with its value 
on the sphere and compatible with the idea that the fractal geometry of 
the disk far from the boundary is identical to the quantum geometry of the 
sphere and the expectation that the influence of the boundary decreases 
as one looks at spins further and further into the disk. This leads to the 
prediction that, in the continuum limit, the bulk quantities should exhibit 
KPZ behavior. 
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Figure 14: Number of vertices at geodesic distance r from the boundary at 
the tricritical point with H = 0.433. 
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Figure 15: Magnetization of vertices at geodesic distance r from the bound- 
ary at the tricritical point with H = 0.433. 
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Figure 16: Sum of spins of vertices at geodesic distance r from the boundary 
at the tricritical point with H = 0.433. 



Specifically, we have examined the scaling of the bulk magnetization and 
susceptibility at the tricritical point. Fig. llTl shows a plot of the logarithm of 
the susceptibility versus the logarithm of the area for a variety of boundary 
magnetic fields. The lines show good power law fits at a variety of H with 
exponents ranging from 0.6 to 0.7 which is in approximate agreement with 
the KPZ prediction of 2/3. 

The situation for the bulk magnetization is more complicated. We would 
expect that this would decrease with area according to the simple power law 
(M^) ~ with A = g predicted by KPZ and consistent with Eq. l47l using 
/ = \/]4. However, the finite size scaling result above as well as fits of the 
same data as used in Fig. ^Jto powers of area yield exponents closer to 0.1. 
This value differs quite markedly from the prediction. This may refiect a 
finite size effect involving the infiuence of spins close to the boundary where 
cIh ~ 2. It is interesting to observe that the resulting exponent for M(r; A'^), 
Adn ~ 0.37 is near the value obtained in |^ IH] of 0.33 using A = g and 
dn = 2. From Eqs. HTIandUHl 

{M^)A^ = [a{H)r'^. (55) 
We hypothesize that the finite size effects can be taken into account by 
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Figure 17: Scaling of bulk magnetic susceptibility with area for > at 
the tricritical point. 



generalizing this to 

l,M^)A'' = f[HA''). (56) 

(Note that rescaling K with A is equivalent to rescaling it with I since we 
have / = \J\A))) Assuming a < and expanding for large A, 

(M^) « hA-"" + fiHA-''+- + ^f2H^A-'+^- + ... (57) 

Since we know that (M^) = for zero H we must have /o = 0. We found our 
data fit this well with a = -0.067 ± 0.005 and & = 0.1 ± 0.02 (See Fig. ITsll . 
If our hypothesis is corect, then our data is consistent with 

as predicted. 

To examine the scaling of the bulk magnetization further, we have cal- 
culated 

"^T^r^ (59) 
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Figure 18: (M^) A^-^ vs. HA'^-^^'^ at the tricritical point. 



which defines an efi'ective (distance dependent) power law exponent for this 
correlation function. Fig. ^Jshows these exponents for r = 1 to 5, averaged 
over several values of H from 0.24 to 1.3. The values for the different areas 
coincide for d = 1 and then decrease at different rates as r increases. The 
rate of decrease is smaller for larger areas lending support to the hypothesis 
that a simple power law decay may indeed describe the behavior in the 
thermodynamic limit. The value of such a power law's exponent can be 
read off as approximately —0.3. This is consistent with the finite size scaling 
results above. 



7 Conclusion 

We have mapped the phase diagram corresponding to a model of Ising spins 
on a fluctuating geometry with disk topology. At zero boundary magnetic 
field H we find three phases separated by critical lines on which either 
the magnetic or boundary length susceptibilities are divergent with disk 
area. Only one of these phase boundary lines corresponds to continuous 
transitions — the line separating ferromagnetic and paramagnetic phases 
on geometries with small mean boundary length. The measured exponent 
for this transition compares well with predictions from Liouville theory. The 
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Figure 19: Exponents of bulk magnetization decaying with distance from 
the boundary at the tricritical point. 



third phase consists of elongated geometries whose boundary length scales 
linearly with disk area. The continuum limit is thought to reside at the 
junction of the these two critical lines — a tricritical point at which the 
spins display long range correlations and the mean boundary size scales as 
the square root of the area (/) ~ \/]4. 

We have tested the results of our simulations in two ways; by comparison 
to hand counted amplitudes for small disks including the Ising sector and 
by a direct computation of the critical exponents governing the pure gravity 
disk partition function. The latter constitutes the first numerical check on 
the values of these exponents and was accomplished by a measurement of the 
baby universe distribution (defined now for disk topology) and the density of 
3-fold coordinated vertices. In both cases good agreement with predictions 
was obtained. 

We have also computed the dependence of the bulk and boundary mag- 
netizations and susceptibilities on the boundary field H at the tricritical 
point. We find the boundary magnetization is independent of boundary 
length or bulk disk area and varies smoothly with H. The corresponding 
susceptibility approaches a constant as H ^ for large area. Both of these 
conclusions are consistent with the predictions of a matrix model calculation 
of a similar model by Carroll, Ortiz and Taylor |3 El E] and at least in the 
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case of the susceptibility differ qualitatively from the situation in flat space. 

In the thermodynamic limit and at the tricritical point, a typical bulk 
spin will be situated infinitely far from the boundary and thus we might 
expect the critical behavior of the bulk magnetization and susceptibility to 
be same as on a sphere independent of the boundary field. Our numerical 
data for the bulk susceptibility is indeed consistent with this — for a wide 
range of fields H we see a power law dependence on the area with exponent 
in the range 0.6 — 0.7 (the continuum/matrix model prediction is |). The 
analysis of the magnetization data is more complicated — a naive power law 
fit of the area dependence at fixed H yields an exponent of approximately 
0.1 which appears quite different from the predicted value of ^. To explain 
this we must assume that the influence of the boundary spins induces large 
finite size corrections in this exponent. 

One conjecture made in the paper of Carroll et al. appears to be false 
— when we measure the bulk magnetization we find that it increases with 
applied field — contrary to the rather surprising result of the matrix model 
calculation where a decrease of bulk magnetization with increasing bound- 
ary field was found. In that case it was conjectured that the physical inter- 
pretation of that result was due to the appearance of a 'neck' in the disk 
geometry which allowed the average bulk spin to be pushed further out from 
the boundary with increasing field. In our case we have verified that the 
opposite is true — the average bulk spin is drawn closer to the boundary 
with increasing field. The explanation for this discrepancy appears to stem 
from a lack of universality in the calculations — our lattice ensemble dif- 
fers from the one used in the analytical calculation and renders detailed 
comparison between the two problematical. Specifically, the function a{H) 
which induces this result is clearly not universal in character. Further tests 
of universality obtained by placing the Ising spins on triangles are currently 
underway. 

On a more positive note we have checked another conjecture made in 

the matrix model paper which speculates that the decay of the boundary 

magnetization into its value in the bulk is governed by a universal power 

of the geodesic distance to the boundary where the power law exponent is 

related to the dressed magnetic exponent of 2-d quantum gravity. We have 

measured both the number of vertices and the magnetization as a function 

of geodesic distance from the boundary and found that they exhibit finite 

size scaling forms for geometries close to the tricritical point. Just as on 

the sphere it appears that these functions reveal the existence of a single 
1 

length scale ^ ~ ^ <^ff governing the behavior of the average number of points 
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and magnetization as we move from the boundary into the interior of the 
disk. Indeed, the fractal dimension dn appears to be close to its value on 
the sphere for critical Ising spins coupled to gravity and lends support to 
the idea that the dressing of Ising spin operators will just follow the KPZ 
predictions. However, we find the extracted value of the bulk magnetic 
operator differs from its exact value which we attribute to large residual 
finite size corrections. 

It would also be interesting to see which of these qualitative features 
carried over to the case of the 3-state Potts model coupled to 2-d gravity on 
the disk. 
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